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ABSTRACT 


The  blade-element  performance  of  a small  transonic 
compressor  rotor  operating  at  high  subsonic  speeds  is 
reported.  Measurements  were  made  downstream  using  a small 
calibrated  pneumatic  probe,  with  simultaneous  measurements 
of  total  pressure,  velocity,  turning  angle,  pitch  angle, 
and  temperature.  In  addition,  a program  is  presented  to 
calculate  blade-element  performance  based  on  empirical  data 
contained  in  NASA  SP-36. 
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The  study  reported  here  is  part  of  an  ongoing  program 
at  the  Naval  Postgraduate  School  to  determine  the  perform- 
ance of  a transonic  compressor  stage.  The  overall  goal  of 
the  program  is  to  develop  the  knowledge  necessary  to  design 
compressors  with  multiple  transonic  stages  which  are  effi- 
cient, light  in  weight,  and  free  of  flutter  problems.  It 
has  long  been  appreciated  that  increasing  the  wheel  speed 
and  relative  air  velocities  in  axial  compressors  would  re- 
sult in  fewer  stages  for  a given  energy  addition.  However, 
in  order  to  achieve  high  overall  efficiencies  in  multiple 
axial  stages,  the  ability  to  predict  the  behaviour  of  a 
single  stage  correctly  must  be  shown.  The  purpose  of  this 
study  was  to  measure  the  blade-element  performance  of  a 
transonic  compressor  rotor  operating  at  high  subsonic  speeds 
and  to  compare  the  measurements  with  the  performance  calcu- 
lated according  to  the  method  given  in  Reference  1. 

The  transonic  compressor  and  drive  turbine  were  de- 
signed by  Dr.  M.  H.  Vavra.  The  program  is  supported  by  the 
Naval  Air  System  Command,  Code  3!0»  through  the  office  of 
Dr.  H.  J.  Mueller. 

To  determine  the  blade-element  performance  of  a com- 
pressor blade  row,  it  is  necessary  to  evaluate  radial  veloc- 
ity and  angle  distributions  into  and  out  of  the  blade  row. 

In  a transonic  compressor,  at  or  near  design  speed,  the 
unsteady  nature  of  the  flows  with  large  pressure  amplitudes 
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complicates  the  problem.  In  addition,  the  possibility  that 
wakes  might  induce  vibrations  in  the  highly  stressed  blades 
prohibit  the  use  of  probes  upstream  of  the  rotor.  Measure- 
ments can  be  made  downstream  of  the  rotor;  however,  since 
the  flow  area  is  small,  boundary  effects  on  probe  measure- 
ments can  be  significant. 

In  order  to  carry  out  the  present  study,  a small  com- 
bination temperature -pneumatic  probe  was  designed  by  Dodge 
(Ref.  2)  specifically  to  be  used  immediately  downstream  of 
the  rotor.  In  previous  work,  the  flow  upstream  of  the  rotor 
was  determined  at  different  flow  rates  (Ref.  3).  It  was 
left  to  the  present  work  to  represent  the  inlet  flow  field 
analytically,  to  calibrate  the  combination  probe  for  bound- 
ary effects,  and  both  to  measure  and  to  calculate  the  rotor 
blade  element  performance. 

In  this  report  the  transonic  compressor  and  test  rig 
are  first  described  in  Section  II.  The  method  of  approach 
and  program  of  measurements  are  described  in  Section  III. 

The  results  of  the  measurements  are  then  described  and 
discussed  in  Section  IV,  and  conclusions  are  given  in  Sec- 
tion V.  In  order  to  maintain  a clear  presentation  of  the 
relevant  results,  the  details  of  the  calculations  and  of 
the  calibration  measurements  are  given  in  separate  Appendices. 
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THE  TRAN 
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COMPRESSOR 


A.  DESCRIPTION 

A drawing  of  the  compressor  test  rig  is  shown  in  Pig.  1. 
The  test  rig  consists  of  an  air-turbine  drive  unit  and  an 
induction  section  containing  a filter,  a throttle,  a set- 
tling chamber  and  a flow  measuring  nozzle.  The  air  turbine 
is  designed  to  deliver  450  horsepower  at  30,000  RPM.  The 
relative  blade  tip  Mach  number  of  the  rotor  at  design  point 
is  1.5*  The  compressor  flow  rate  is  controlled  by  an  elec- 
tric-hydraulic rotating  throttle  plate.  Power  is  provided 
by  a compressed  air  supply  from  an  Allis-Chalmers  multi- 
stage axial  compressor.  A complete  description  of  the 
Turbopropulsion  Laboratory  and  test  facilities  is  contained 
in  Ref,  4.  A drawing  of  the  transonic  compressor  is  shown 
in  Fig.  2. 

The  blade  design  for  the  transonic  compressor  rotor  was 
derived  from  circular-arc  elements  wrapped  on  six  conical 
surfaces.  The  outer  surface  (at  the  case  wall)  is  a cylin- 
der. The  blading  geometry  is  given  in  Fig.  3.  The  blade 
pressure  side  is  flat  and  the  suction  side  is  a circular- 
arc.  The  solidity  is  fairly  constant  with  a value  of  1.35 
at  the  mid-streamline.  The  outer  diameter  of  the  rotor 
blade  is  11  inches  with  a channel  height  of  1.92  inches 
downstream  and  a radius  ratio  at  the  inlet  face  of  0.5. 
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B.  INSTRUMENTATION 

The  transonic  compressor  is  equipped  with  stationary 
probes  to  provide  routine  performance  data.  Triple  redun- 
dancy is  provided  in  the  measurement  of  the  power j either 
the  mass  flow  rate  and  temperature  drop  through  the  drive 
turbine,  or  the  mass  flow  rate  and  temperature  increase 
through  the  compressor,  or  the  speed  and  the  measurement  of 
torque  can  be  used.  The  torque  is  measured  by  strain  gauges 
on  flexures  against  which  the  stator  assembly  is  free  to 
rotate.  Flow  straighteners  before  and  after  the  stage  en- 
sure that  the  stator  and  rotor  torques  are  equal.  The 
method  based  on  torque  and  speed  is  the  most  accurate  since 
it  is  independent  of  flow  profiles.  Static  pressure  taps 
were  installed  on  both  outer  and  inner  walls  of  the  com- 
pressor casing.  The  rotative  speed  of  the  test  rotor  was 
determined  by  a magnetic  flux  cutter  and  electronic  counter. 
The  fixed  instrumentation  will  not  be  described  further. 

It  was  used  in  the  present  tests  only  to  provide  the  refer- 
ence conditions  for  the  survey  measurements.  The  compressor 
cylindrical  casing  was  modified  to  accommodate  a variety 
of  probes  and  sensors  which  may  be  required  in  the  program 
of  steady  state  and  unsteady  flow  pressure  measurements. 

The  casing  can  be  rotated  to  position  probes  in  the  periph- 
eral direction. 


The  rotor  exit  survey  measurements  reported  here  were 
obtained  using  a combination  probe  (Fig.  4)  designed  by 
Lt.  F.  J,  Dodge  (Ref.  2).  The  combination  probe,  following 
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calibration,  was  used  to  determine  total  pressure  rise,  to- 
tal temperature  rise,  and  flow  velocity  in  magnitude  and 
direction.  In  addition,  static  pressure  was  measured  using 
the  wedge  static  probe  Shown  in  Fig.  5»  which  was  based  on 
a NASA  design  provided  by  Mr.  Lloyd  N.  Krause  of  NASA  Lewis 
Research  Center.  Each  probe  had  two  "static  pressure"  sen- 
sors Whi bh,  if  the  probes  were  rotated  until  the  indicated 
static  pressures  were  equal,  allowed  the  probes  to  be  aligned 
with  the  flow  direction.  It  was  determined  that  better 
accuracy  could  be  attained,  however,  by  aligning  both  probes 
to  the  flow  direction  sensed  by  the  combination  probe.  This 
was  because  the  sensors  on  the  wedge  probe  were  separated 
in  the  radial  direction.  A chrome 1-constantan  fine-wire 
thermocouple  in  the  combination  probe  was  used  in  conjunction 
with  a similar  sensor  mounted  upstream  of  the  rotor  to 
measure,  differentially,  the  total  temperature  rise. 

Radial  surveys  were  taken  at  a fixed  peripheral  case 
angle  by  manually  positioning  the  probe  in  the  flow  using  a 
United  Sensor  traverse  unit,  manually  rotating  the  combina- 
tion probe  to  balance  indicated  static  pressures,  and  ad- 
justing the  wedge  probe  to  the  same  angle.  All  pressures 
were  measured  by  a calibrated  transducer  in  a Scanivalve 
arrangement.  (A  capability  exists  to  remotely  control  probe 
movement  and  balancing  through  the  use  of  a programmable 
microprocessor  and  a probe  traversing  mechanism  designed  by 
Lt,  D,  D,  Patton  (Ref,  5).  Remote  operations  will  be  needed 
as  higher  operational  speeds  are  developed.) 
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C.  DATA  ACQUISITION  AND  REDUCTION 

Pressures,  temperatures,  rotational  speed  and  other  data 
were  scanned  and  punched  onto  paper  tape  by  a B & F Model 
24133  solid  state  data  acquisition  system.  A complete  de- 
scription of  this  system  is  contained  in  Ref.  4.  The 
punched  tape  data  was  processed  using  Hewlett-Packard  Model 
983OA  programmable  calculator  data  system.  The  data  was 
transferred  to  a Model  9867B  mass  memory  disk  unit  for  stor- 
age prior  to  data  reduction.  The  data  facilities  are  shown 
in  Fig.  6. 

A program,  designated  "WH001"  was  written  to  transfer 
raw  data  from  tapes  to  the  mass  memory.  A second  program 
"HAWK1"  was  written  to  calculate  temperature,  corrected 
non-dimensional  velocity,  turning  angle,  total  pressure, 
and  static  pressure.  Two  other  programs,  used  in  sequence 
("REPR01"  and  "RA5C" ) completed  the  data  reduction  process. 
Plotting  routines  developed  by  Dr.  R.  P.  Shreeve  were 
adapted  to  plot  distributions  of  Mach  number,  pressure  rise, 
temperature  rise,  flow  angle,  and  loss  coefficients,  using 
the  X-Y  plotter. 

In  reducing  the  probe  measurements  to  obtain  velocities 
and  loss  data,  the  thermocouple  calibration  reported  in 
Ref.  2 was  used.  However,  new  pneumatic  calibrations  were 
performed  and  corrections  were  derived  as  described  in 
Appendix  B. 


III.  METHOD  OF  APPROACH 

I 

A.  METHOD 

The  usual  approach  to  the  design  of  axial  compressors  is 
described  in  Ref.  1.  The  method  involves  combining  an  axi- 
symmetric  through-flow  analysis,  which  prescribes  the  stream 
surfaces , with  a two-dimensional  cascade  performance  calcu- 
lation which  determines  the  required  blading  geometry  on 
those  surfaces.  The  blading  is  generally  layed  out  on  coni- 
cal approximations  of  the  prescribed  stream  surfaces. 

When  a new  machine  is  built,  the  opportunity  is  there 
to  compare  the  actual  blade  element  performance  established 
by  measurement  with  that  which  can  be  calculated  by  applying 
the  design  approach.  The  comparison  can  be  used  to  derive 
corrections  for  use  in  future  designs,  and  it  is  particularly 
useful  to  obtain  data  for  new  blade  profiles. 

The  transonic  compressor  investigated  here  was  designed 
using  the  above  approach,  and  the  resulting  blading  is  shown 
in  Pig.  3.  It  is  to  be  noted  that  the  rotor  blading  is  not 
typical  in  design.  In  order  to  establish  the  blade  element 
performance  of  a rotor  experimentally,  measurements  should  • 
be  made  of  the  flow  into  and  out  of  the  rotor  on  the  stream 
surfaces  on  which  the  blading  profiles  are  prescribed.  How- 
ever, it  is  desirable  not  to  operate  at  high  wheel  speeds 
with  strong  wakes  from  upstream  probes,  and  so  in  the  pre- 
sent work  a method  of  prescribing  the  inlet  flow  field  from 
measurements  of  only  total  flow  rate  and  wheel  speed  was 
derived.  The  method  is  given  in  the  following  section. 


The  flow  downstream  of  the  rotor  was  determined  using  two 
calibrated  probes.  The  calibration  and  corrections  re- 
quired for  the  probes  are  described  in  Appendix  B.  Prom 
the  rotor  exit  flow  measurements,  the  blade  element  loss 
coefficient  and  the  local  deviation  angle  were  calculated. 

Finally,  as  described  in  Appendix  A.  the  method  of 
calculating  deviation  angle  and  loss  coefficient  given  in 
Chapter  7 of  Reference  1 was  programmed  for  the  Hewlett- 
Packard  Model  983OA  calculator  so  that  these  quantities 
could  be  computed  for  four  blade  element  surfaces  at  any 
given  compressor  flow  rate  and  operating  speed. 

B.*  INLET  FLOW  FIELD  REPRESENTATION 

The  flow  field  ahead  of  the  rotor  was  measured  in  early 
tests  using  a calibrated  pneumatic  probe  (Ref.  3).  The 
measurements  were  compared  with  the  results  of  a compress- 
ible finite  difference  calculation  reported  in  Reference  6, 
and  blockage  factors  were  derived  to  relate  the  inviscid 
calculations  to  the  measured  flow  field. 

The  method  adopted  here  was  based  only  on  Anderson's 
measured  data.  As  described  in  Appendix  C a polynomial  re- 
lation was  found  to  correlate  the  measured  velocity  distri- 
butions away  from  the  hub  and  tip  regions.  The  flow  function 
in  this  region  was  therefore  calculated  using  Eq.  C(3)  and 
Eq.  C(7) . 
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C.  USE  OF  FLOW  PROBES 

Pressure  probes  are  designed  to  exploit  the  distribution 
of  pressure  which  occurs  over  a body  immersed  in  a moving 
fluid.  These  pressure  variations  depend  mainly  on  the  flow 
velocity  so  that  with  a carefully  selected  choice  of  body 
shape,  having  multiple  sensors^a  probe  can  be  calibrated  to 
determine  flow  velocity.  The  relationship  between  pressure 
and  velocity  must  be  established  over  a range  of  Mach  num- 
bers. In  addition  to  flow  velocity  in  magnitude  and  direc- 
tion, the  distributions  of  stagnation  temperature,  and  total 
and  static  pressure  are  usually  required.  The  probe  designed 
by  Dodge  and  reported  in  Ref.  2 is  a small  lightweight  probe 
of  aerodynamic  design  which  allows  these  quantities  to  be 
determined  at  any  position  in  the  flow. 

In  order  to  provide  an  independent  measurement  to  com- 
pare with  the  combination  probe,  the  static  pressure  distri- 
bution was  measured  with  a second  probe,  and  static  pressures 
were  also  recorded  at  hub  and  tip  wall  taps.  The  wedge 
static  probe  shown  in  Fig.  5 has  the  favorable  characteristic 
that  when  it  is  aligned  with  the  flow  in  the  yaw  plane  (by 
balancing  the  pressures  measured  on  the  wedge  surfaces)  the 
indicated  static  pressure  is  relatively  insensitive  to  angle 
variation  in  the  pitch  plane.  As  reported  in  Ref,  7 this 
type  of  probe  has  been  used  successfully  up  to  a Mach  number 
of  1.6.  Further,  for  measurements  in  the  transonic  regime 
the  probe  has  favorable  characteristics.  If  the  leading 
edge  is  sharp  and  the  wedge  is  small,  Ref.  7 states  that  the 
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calibration  is  also  independent  of  Mach  number.  This  was 
found  to  be  the  case  for  the  probe  used  in  the  present  work. 

D.  PROGRAM  OF  MEASUREMENTS 

The  combination  probe  and  the  wedge  static  probe  were 
mounted  at  the  exit  of  the  rotor  with  a peripheral  separa- 
tion of  120  degrees.  Four  surveys  were  conducted  at  50 
percent  of  design  speed,  at  throttle  conditions  varying 
from  near  surge  to  maximum  flow  rate.  Case  angle  for  the 
first  survey  was  3^5  degrees.  Subsequent  surveys  were  con- 
ducted at  a case  angle  of  357  degrees  because  it  was  deter- 
mined that  peripheral  variations  induced  by  the  stator 
blading  were  felt  least  at  this  setting.  RPM  was  the  same 
(15,230130)  for  all  four  runs. 

The  procedure  in  each  survey  was  to  move  the  probes  in- 
ward simultaneously,  balance  the  combination  probe  to  sense 
the  flow  angle,  and  then  set  the  wedge  probe  to  the  measured 
angle.  A data  scan  was  taken,  and  the  procedure  was  repeated 
at  0.1  inch  intervals.  Because  of  the  hooked  geometry  of 
the  wedge  probe,  measurements  with  this  probe  of  static 
pressure  could  not  be  made  closer  than  .6  inches  to  the 
outer  wall  boundary.  Static  pressures  were  derived  for  the 
complete  passage  through  the  use  of  the  combination  probe 
calibration.  Wall  static  pressure  measurements  at  the  hub 
and  tip  were  also  recorded.  Yaw  angle  and  linear  displace- 
ment were  recorded  manually.  All  other  data  were  automati- 
cally recorded  by  the  data  recording  system.  Pressures  were 
recorded  to  an  accuracy  of  .1  inches  of  water. 
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V.  RESULTS  AND  DISCUSSION 

A.  TIME-AVERAGED  PLOW  FIELD 

The  distributions  of  total  pressure  rise,  total  temper- 
ature rise,  Mach  number  and  flow  angle  measured  using  the 
combination  probe  downstream  of  the  rotor  are  shown  in 
Figs.  7,  8,  9,  and  10. 

The  total  pressure  rise  was  nearly  constant  with  radius, 
whereas  the  temperature  rise  was  indicated  to  be  larger 
toward  the  tip.  The  larger  indicated  temperature  rise  near 
the  case  wall  might  be  caused  by  higher  temperatures  in  the 
incoming  case  wall  boundary  layer.  Both  the  Mach  number 
and  the  flow  angle  increase  toward  the  hub,  so  that  the 
through-flow  velocity  tends  to  remain  constant.  The  results 
for  the  flow  angle  are  qualitatively  consistent!  however,  a 
large  decrease  can  be  observed  in  the  magnitude  of  the  flow 
angle  for  the  open  throttle  case  (Fig.  8).  No  definitive 
explanation  was  found,  but  as  will  be  shown  later  this 
measurement  leads  to  inconsistencies  in  the  correlation  of 
blade  element  performance. 

The  velocity  and  pitch  angle  distributions,  which  were 
derived  from  the  probe  calibration  and  corrected  for  bound- 
ary effects  (Appendix  B),  are  shown  in  Fig.  11  and  Fig.  12. 
The  pitch  angle  shows  a smooth  variation  between  the  values 
of  the  wall  slopes  at  the  hub  and  tip  boundaries.  It  is 
shown  in  Appendix  B that  the  pitch  angle  is  determined  to 
only  l£°  accuracy,  so  that  the  variations  in  pitch  angle 
indicated  in  Fig.  12  might  be  exaggerated. 
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A sensitive  test  of  the  accuracy  of  the  combination  probe 
measurements  was  to  examine  the  static  pressure  distribution 
inferred  from  the  measured  velocity  and  probe  impact  pres- 
sure measurement.  Pig.  13  shows  the  distribution  of  static 
pressure  calculated  in  this  way.  Also  included  in  Fig.  13 
are  the  values  of  static  pressure  measured  at  the  hub  and 
tip  wall  static  taps.  It  can  be  seen  that  agreement  near 
the  case  wall  was  excellent.  The  comparison  near  the  hub 
wall  was  not  as  good?  however,  very  similar  results  were 
obtained  from  the  wedge  probe  as  shown  in  Fig.  14.  While 
no  readings  were  possible  near  the  outer  wall  using  the 
wedge  probe,  the  two  probes  gave  very  similar  results  near 
the  hub  wall,  and  gave  almost  identical  results  for  the  sur- 
vey at  open  throttle. 

As  shown  in  Appendix  3,  the  effect  of  flow  boundaries  on 
the  readings  of  the  combination  probe  is  small.  The  magni- 
tude of  the  correction  which  was  included  in  the  above  re- 
sults can  be  seen  in  Fig.  15.  The  correction  is  less  than 
1 % in  velocity  over  the  useful  range  of  the  survey. 

Velocity  diagrams  near  the  mean  radius  (0.55  of  the 
channel  height)  are  shown  in  the  four  parts  of  Fig.  16. 

B.  DEVIATION  ANGLES 

The  design  procedure  for  calculating  the  deviation 
angles  for  a given  blade  profile,  with  corrections  for 
three-dimensional  effects,  is  described  in  Appendix  A. 

Blade  profiles  are  known  on  four  conical  surfaces  in  the 


flow  field,  on  which  the  blading  dimensions  are  those  given 


in  Fig,  3.  Using  relative  flow  angle  from  the  velocity 
diagram  data  (Fig,  16  for  example),  the  flow  deviation  angle 
(S)  can  be  calculated  from  the  stagger  angle  (TO  and  camber 
angle  (£)  using  the  relation 


The  results  for  the  four  probe  surveys  are  shown  in 
Fig,  17 -in  comparison  with  the  results  of  calculations 
following  Ref.  1. 

Departures  from  values  predicted  using  Ref.  1 are  quite 
large.  However,  the  design  procedure  strictly  applies  only 
to  the  case  where  the  operation  is  at  the  point  of  minimum 
loss  for  the  given  blade  element,  A closer  correlation 
might  be  expected  as  the  higher  operating  levels  are  achieved 
and  closer  to  design  conditions  are  attained,  (Ref.  1 states 
that  at  transonic  inlet  Mach  number  levels,  the  point  of 
minimum  loss  should  constitute  the  desired  design  setting, ) 

It  is  noted  that  the  through-flow  velocity  is  considerably 
increased  through  the  rotor  at  the  50%  design  speed  setting, 
and  this  would  effect  the  deviation  angle.  No  explanation 
was  found  for  the  indication  of  negative  values  of  deviation 
angle  near  the  hub,  although  a careful  check  was  made  on  the 
accuracy  of  the  probe  yaw  angle  measurements  and  on  the  ac- 
curacy of  the  data  reduction  program.  It  must  tentatively 
be  concluded  that  there  are  strong  three-dimensional  effects 
present  which  can  lead  to  the  indication  (locally)  of  nega- 
tive deviation  angles.  More  data  must  be  obtained  to  investi- 
gate this  speculation. 


C . LOSSES 

Loss  coefficients  for  the  flow  over  the  blade  elements, 
assuming  that  the  stream  surfaces  were  the  conical  surfaces 
given  in  Fig.  3»  were  calculated  from  the  probe  measurements 
using  the  methods  described  in  Appendix  D.  The  calculated 
loss  coefficient  (&')  was  consistent  with  the  definition 
of  loss  coefficient  used  in  Ref.  1.  Of  the  two  methods 
described  in  Appendix  D,  Method  1 uses  both  the  stagnation 
pressure  rise  and  stagnation  temperature  rise  measurements 
of  the  probe.  Method  2 does  not  use  the  stagnation  temper- 
ature rise  measurement,  but  uses  a temperature  rise  calcu- 
lated from  the  velocity  and  angle  measurements,  through 
Euler's  turbine  equation. 

The  results  are  shown  in  Fig.  18.  Good  agreement  was 
obtained  between  the  two  methods  of  calculating  the  loss 
coefficient,  however  negative  values  were  calculated  for  the 
inner  stream  surface.  It  should  be  noted  that  all  measured 
data  and  all  blade  properties  varied  smoothly  with  radius, 
so  that  smooth  Gurves  between  the  individual  data  values 
could  be  drawn.  The  distributions  are  in  qualitative  agree- 
ment with  earlier  measurements,  and  it  is  noted  that  nega- 
tive values  of  losses  were  previously  reported  by  Dodge 
(Ref.  2)  from  measurements  made  with  the  same  combination 
probe. 

A negative  value  of  the  loss  coefficient  can  result  only 
from  a measurement  of  stagnation  pressure  rise  which  is  too 
high,  or  a measurement  of  stagnation  temperature  rise  which 
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is  too  low.  This  can  he  seen  if  Eq.  D(13)  is  rewritten  in 
the  form 
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wherelj^is  the  ratio  of  total  temperatures,  (T^g/Tt^). 
Clearly,  since  XUg  and  XU1  are  positive, <Jut  is  negative  only 
if  > (Pt2/pti^ • Since  the  prediction  of  Method  1 and 


Method  2 agree,  and  because  a favorable  comparison  with  ear- 
lier unreported  measurements  of  temperature  rise  was  made, 
it  is  tentatively  concluded  that  the  impact  pressure  rise 
measurement  is  too  large.  A positive  error  in  impact  pres- 
sure could  result  from  the  pneumatic  integration  of  the  peri 
odic  unsteadiness  of  the  flow  from  the  rotor.  All  parts  of 
the  probe  are  the  open  ends  of  continuous  tubing.  It  is 
known  that  pneumatic  probes  of  different  designs  integrate 
differently  in  periodic  flow  conditions.  An  investigation 
of  the  pneumatic  averaging  characteristics  of  the  probe  is 
required  in  order  to  correct  the  measurements  of  losses. 

It  should  be  noted  that  the  velocity  and  pitch  angle  calcu- 
lated from  the  time -averaged  probe  readings  should  not  be 
as  sensitive  to  the  averaging  behaviour  of  the  probe.  This 
is  because  all  parts  are  made  from  similar  tubing  and  will 
average  in  a similar  way.  The  loss  coefficient,  however. 


is  very  sensitive  to  errors  in  impact  pressure  measurement 


While  the  magnitude  of  the  loss  coefficients  is  clearly 
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in  error,  it  is  of  interest  to  attempt  to  correlate  the 
measurements  on  the  basis  of  diffusion  factor  as  defined 
Ref.  1.  The  result  of  this  attempt  is  shown  in  Fig.  19. 
diffusion  factor  (D)  was  calculated  from  the  definition 


in 

The 


D = 1 


W2  + R2WU2 

Wl  O' (R1  + R2) 


where  W denotes  the  relative  velocity,  R is  the  radius,  (5"  is 
the  solidity  and  subscripts  1 and  2 denote  stations  upstream 
and  downstream  of  the  rotor.  Wu  is  the  peripheral  component 
of  the  relative  velocity.  The  "total  loss  parameter", 

( co^cos  0 2/2<T ) is  consistent  with  that  shown  in  Fig.  203 
of  Ref.  1 and  was  calculated  from  probe  measurements. 

It  can  be  seen  in  Fig.  19  that  a reasonable  correlation 
of  the  loss  coefficient  data  on  the  basis  of  diffusion  factor 
was  obtained  for  three  of  the  four  sets  of  data.  The  data 
for  full  open  throttle  departed  significantly  from  the  data 
from  three  other  test  conditions.  It  should  now  be  recalled 
that  the  flow  angles  measured  for  this  test  were  also  incon- 
sistent. More  data  is  again  needed,  and  a correction  must 
be  devised  for  unsteady  flow  effects  on  the  probe  measure- 
ments before  a complete  correlation  can  be  obtained. 
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VI.  CONCLUSIONS 
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The  attempt  was  made  here  to  determine  the  blade  element 
performance  of  an  "impulse"  rotor  designed  for  transonic  flow 
(design  point)  but  operating  at  high  subsonic  (off -design) 
conditions.  The  results  were  preliminary,  but  served  to 
define  the  problems  involved  in  obtaining  necessary  design 
data  from  empirical  correlations  of  measurements  made  in 
actual  machines. 

The  results  obtained  depended  on  the  performance  of  a 
particular  combination  probe  when  used  downstream  of  the 
ro'cor.  The  following  conclusions  were  drawn  concerning  the 
probe  measurements  i 

(i)  The  effect  of  the  presence  of  flow  boundaries  on  the 
probe  measurements  was  small,  but  corrections  were  derived 
and  applied  in  the  present  work.  Static  pressure  derived 
from  the  probe  measurements,  static  pressure  from  wedge  probe 
measurements  and  wall  static  pressures  were  in  reasonable 
agreement. 

(ii)  The  temperature  rise  across  the  rotor  measured  by 
the  probe,  using  the  probe  calibration,  was  consistent  with 
the  rise  deduced  from  Euler's  turbine  equation.  All  indica- 
tions were  that  the  probe  provided  an  accurate  measurement 
of  the  temperature  rise. 

(iii)  The  impact  pressure  measured  by  the  probe  in  the 
unsteady  flow  downstream  of  the  rotor  was  probably  too  high. 
An  investigation  of  the  unsteady  response  of  the  tubing  used 
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in  the  probe  design  is  required  in  order  to  correct  the 
data  obtained  so  far. 

(iv)  The  flow  velocity  derived  from  the  probe  measure- 
ments was  probably  not  as  sensitive  to  unsteady  flow  effects 
since  all  the  pressure  sensors  on  the  probe  were  constructed 
from  the  same  tubing.  The  derived  velocity  depends  mainly 
on  differences  between  the  pressures  registered  by  the  sen- 
sors. In  steady  flow  the  velocity  was  determined  to  better 
than  1.0#  accuracy,  to  nearly  sonic  speeds. 

(v)  The  pitch  angle  was  determined  by  the  same  probe  to 
within  l£°. 

(vi)  The  wedge  probe  provided  a good  measurement  of 
static  pressure  in  the  compressor  when  rotated  to  a pre- 
measured flow  angle. 

As  a result  of  (iii),  the  loss  coefficients  determined 
from  probe  measurements  were  too  small,  becoming  negative 
over  the  blading  near  to  the  hub.  A correction  to  the  im- 
pact pressure  measurements  and  more  data  are  needed  before 
a useful  correlation  of  the  losses  can  be  obtained.  However, 
it  was  found  that  the  present  results  from  three  tests  and 
four  different  blade-element  surfaces  correlated  reasonably 
well  in  terms  of  the  NASA  diffusion  factor. 

The  deviation  angles  measured  in  the  present  tests  de- 
parted significantly  from  the  values  predicted  using  NASA 
SP-36.  No  explanation  was  found  for  the  negative  values 
measured  near  to  the  hub,  although  the  possibility  exists 
that  they  could  be  caused  by  three-dimensional  effects  in  the 
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flow  through  the  rotor.  More  data  are  needed  to  confirm  the 
trends  in  the  deviation  angles.  Agreement  with  the  predicted 
values  was  not  expected  since  the  blade  profiles  were  unusual, 
and  were  operating  well  away  from  the  design  condition. 
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PIG.  6 SCHEMATIC  OP  THE  HEWLETT-PACKARD  MODEL  9830A  CALCULATOR  DATA  REDUCTION  SYSTEM 


USING  THE  COMBINATION  PROBE 


NDRMHUZED  DlSfUKEKNT  (HUB  TO  TIP) 

PIG.  13  STATIC  PRESSURE  DISTRIBUTIONS  DOWNSTREAM  OP  THE  ROTOR 
DERIVED  USING  THE  COMBINATION  PROBE  CALIBRATION  AND 
IMPACT  PRESSURE  MEASUREMENT 


NDRHHLiZED  DI5PUKEMCNT  (NUB  TO  TIP) 

FIG.  14  STATIC  PRESSURE  DISTRIBUTIONS  MEASURED  DOWNSTREAM  OF 
THE  ROTOR  USING  THE  WEDGE  STATIC  PROBE 


« PROBE  DIMETERS  FROM  OUTER  HU  BDUNDRRY  OF  THOMSON  1C  COMPRESSOR 

FIG.  15  BOUNDARY  EFFECT  CORRECTION  APPLIED  TO  THE  COMBINATION  PROBE  RESULTS 


6a  VELOCITY  TRIANGLES  AT  MIDSTREAM  FOR  RUN  50  POINT  3 


PIG.  l6b  VELOCITY  TRIANGLES  AT  MIDSTREAM  FOR  RUN  52  POINT  3 


VELOCITY  TRIANGLES  AT  MIDSTREAM  FOR  RUN  52  POINT  23 


FIG.  l6d  VELOCITY  TRIANGLES  AT  MIDSTREAM  FOR  RUN  52  POINT  46 


NDRMnL  I XtCl>  t>  I HPUHCCMCNT  < HUM  TO  TIP) 

FIG.  17  DEVIATION  ANGLES  FROM  PROBE  MEASUREMENTS  AND  CALCULATIONS 


oirruBiaN  rwcTDR 

FIG.  19  DIFFUSION  FACTOR  CORRELATION  OF  LOSS  COEFFICIENTS 


Flow  through  the  blading  of  an  axial-flow  compressor  is 
an  extremely  complicated  three-dimensional  process.  As  a 
result,  simplified  design  models  have  been  devised  that  per- 
mit the  use  of  two-dimensional  techniques.  The  approach  is 
described  fully  in  Reference  1.  Involved  is  first  the  cal- 
culation of  the  "through-flow"  assuming  axial  symmetry,  then 
the  calculation  of  the  "blade-element  flow"  assuming  approx- 
imately two-dimensional  behaviour.  The  latter  calculation, 
based  on  Chapters  6 and  7 of  Reference  1,  is  the  subject  of 
this  section. 

Five  quantities  shown  in  Figure  A1  determine  the  aero- 
dynamic behaviour  of  a cascade  in  ideal  flow.  Three  of  these 
define  the  geometry  of  the  blade » they  are,  the  shape  of  the 
blade,  the  orientation  of  the  blades  with  respect  to  the  cas- 
cade axis,  and  the  solidity.  The  fourth  quantity  identifies 
the  direction  of  flow  ahead  of  the  cascade  and  finally,  the 
Mach  number  must  be  known.  When  these  quantities  are  speci- 
fied, the  inviscid  flow  may  be  determined. 

Real  cascade  flows  however,  experience  many  other  in- 
fluences including  friction,  turbulence,  three-dimensional 
and  unsteady  effects,  and  boundary  layer  interaction.  The 
procedure  outlined  in  Reference  1 correlates  cascade  and  com- 
pressor blade  element  performance  data.  Various  correction 
factors  are  introduced  which  account  empirically  for  the 
effects  which  can  not  be  predicted  analytically.  The  result 
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is  a method  of  calculating  the  performance  of  compressor 
blading.  The  method  can  be  applied  to  both  rotors  and  stators 
over  a fairly  wide  range  of  Mach  numbers.  The  method  is 
limited  however  to  the  calculation  of  performance  close  to 
design  point. 

NASA  SP-36  has  become  a reference  text  in  the  axial 
compressor  industry.  A brief  summary  of  the  method  is  pre- 
sented here,  followed  by  a computerized  solution  which  uses 
the  extensive  correlations  of  data  presented  in  this  manual. 

The  main  thrust  was  to  develop  a program  to  calculate 
the  blade  element  performance  of  a blade  design  of  arbitrary 
geometry.  The  program  would  provide  a very  powerful  design 
tool,  and  would  provide  the  basis  for  a program  to  predict 
the  performance  of  a transonic  compressor  rotor.  Reference  6 
reports  an  extension  of  this  was  used  with  considerable 
success  up  to  a Mach  number  of  about  1.5. 

A2,  METHOD  OF  BLADE  ELEMENT  PREDICTION 

A2.1  Two-Dimensional  Incidence  Angle.  Any  cascade  of 
aerodynamic  blading  behaves  qualitatively  as  shown  in  Fig- 
ure A2  where  A i is  the  "range"  of  performance  and  the  mid- 
point of  the  curve  is  defined  as  the  reference  minimum  loss 
incidence  angle  (iref  )»  For  blading  with  a given  camber 
(0*),  thickness,  distribution,  and  maximum  thickness  ratio 
(t/c),  the  solidity  (<f)  and  blade  angle  ('i ) complete  the 
definition  of  the  cascade  geometry .'  Given  $ ^ and  the  inlet 
velocity  the  method  in  Reference  1 allows  the  determination 
minimum  loss  incidence  angle,  the  deviation  angle,  and  the 
loss  coefficient. 


The  reference  incidence  angle  for  any  blade  element  can 


be  calculated  from  correlations  of  two-dimensional  cascade 
measurements.  Two-dimensional  reference  incidence  angle  for 
1095  thick  blades  is  given  as 

l2.D  * WlO  + **'  A(1) 

where  is  a function  of  blade  thickness  distribution  and 
(iQ)^Q  is  the  zero-camber  incidence  angle  (function  of  ^ ^ 
and  & ),  n is  a slope  factor  given  by 

n = ( i - iQ)  / A(2) 

which  is  also  a function  of^  and  O'  . 

In  Ref.  1,  ( iQ) is  §iven  in  Fig.  1 37  and  n is  given 
in  Pig,  138.  has  the  value  0,7  for  double  circular-arc 
aerofoil  shapes, 

A2.2  Two-Dimensional  Deviation  Angle.  As  in  the  analy- 
sis of  reference  incidence  angle,  the  correlation  of  blade- 
element  deviation  angle  is  expressed  in  terms  of  blade 
geometry  by 

$%-D  * KS  (S0°>10  *im  /<rV'  AO) 

where  K^is  a function  of  blade  type  (0.7  for  double  circular 
arc).  (S0°>io  zero”camber  deviation  angle  (function 

of?i.  and^),  m is  a function  of$lt  for  the  different  basic 
camber  distribution,  and  b is  an  exponent  that  is  also  a 
function  of  ^ 

Since  deviation  angle  varies  somewhat  with  incidence 
angle,  a corrected  deviation  angle  is  defined  as 

^2-D°  = \ (So0)10  + ^c  “ ^-D^^dT^-D 


where  (d£/d^)2_D  is  the  slope  of  the  two-dimensional  varia- 
tion of  deviation  angle  with  incidence  angle  at  reference 
incidence,  and  ift  is  the  actual  incidence  angle. 

Camber  and  turning  angle  are  related  by  the  equation! 


r ’ H -h  + s + 1 


A(5) 


Once  the  two-dimensional  deviation  and  incidence  angles  have 
been  established,  the  turning  angle  can  be  computed. 

A2.3  Corrections  for  Three-Dimensional  Effects.  The  next 
step  in  the  procedure  is  to  correct  reference  deviation  and 
reference  incidence  angles  for  three-dimensional  effects. 

This  is  accomplished  by  defining  the  three-dimensional  refer- 


ence angles 


(iref)3-D  = *2-D  * (i  " i2-D)deS 


A(6) 


(Sref}3-D  * £ 2-D  + (£  " S2-D}  A(7) 

where  (i  - i2_D)  and  ($  - are  functions  of  Mach  number. 

In  summary,  the  desired  turning  angle  is  given  by  Eq. 
288  of  Ref.  1. 

■» 

A ? * -■%♦“]-  $1  " S2-D}  ‘ (ic  ‘ i2-D) 

[l.°  - (^-)2.D]  - Ki^o^lO  + K5^0^10  A(8) 


A(8) 


The  loss  coefficient  for  a blade 


element  can  be  calculated  from  the  correlation  given  in 


Pig.  203  of  Ref.  1,  Here  the  "total  loss  parameter"  is  cor- 
related in  terms  of  "Diffusion  Factor".  The  calculation  of 
losses  is  given  in  Appendix  D. 


A3.  POLYNOMIAL  REPRESENTATION  OP  CORRELATIONS 

All  values  on  the  right  side  of  Eq.  A(8)  may  be  deter- 
mined using  data  given  in  Figures  in  Ref.  1 once  the  geometry 
of  the  blade  has  been  given  and  the  inlet  velocity  calculated. 

Each  of  the  unknown  values  of  Eq.  A(8)  are  functions  of 
blading  geometry  and  air  inlet  angle  or  blading  geometry  and 
Mach  number,  and  are  described  by  well-behaved  families  of 
curves.  Fig.  168  of  Ref.  1 has  been  chosen  to  illustrate 
the  process  of  computerizing  the  solution,  which  involved  re- 
ducing the  data  correlations  given  in  figures  to  analytical 
expressions. 

Eleven  values  of  air  inlet  angle  ( ) were  chosen  for 
each  solidity  given  in  Fig.  168,  The  values  were  plotted 
as  slope  factor  (m)  vs.  air  inlet  angle  ((Jj.)  for  constant 
solidity.  The  data  points  for  each  solidity  were  fit  by  a 
third  order  polynomial  curve  using  the  Hewlett-Packard  9830A 
calculator  and  X-Y  plotter.  For  simplicity,  a fifth  order 
plotting  routine  was  always  used  with  the  higher  order  co- 
efficients equal  to  zero  if  a lower  order  resulted  in  a more 
appropriate  fit.  The  plot  and  polynomial  curve  fits  for 
each  solidity  are  shown  in  Fig.  A3.  The  polynomial  curve 
fit  for  each  set  of  data  points  gave  six  soefficients  for 
each  curve.  Thus  for  each  m vs.  (5  plot,  the  six  coeffi- 
cients AQ,  A1# ,A,-  were  calculated  for  the  polynomial 
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l 

l 


i i 


expression  for  m as  a function  of  ^ ^ i 

m = A0  ♦ A161  + A2012  + A3(Ji3  + A^4  + A^5  A(ll) 

at  each  value  of  solidity  with  converted  to  radians  to 
minimize  the  error  in  higher  order  terms.  The  coefficients 

Aq,  Alf ,A^  derived  from  the  curve  fits  of  Fig.  A3  were 

plotted  against  each  of  twelve  values  of  solidity  as  follows* 


Ao 

(i) 

vs. 

<r  (I).... 

. . . ,12points 

A1 

(i) 

vs. 

<T(I).... 

. . . .12points 

i, 

(i) 

vs. 

(S’  (I).... 

. . . ,12points 

A(12) 

This  produced  curves  of  the  form 

Aq  = aQ  + a^  +3.2^  + a3<r3  + a4<r4  + a^ 

1 1 12  12  1 k 1 < 

A^  = Sq  + 3-^d"  Sg  (S  a^  & a^  a^  (f 


. A(13) 

A5  = aQ5  + a1V  + a25^  + a^3  + 4 + a^3 

which  resulted  in  polynomial  curve  fits  shown  in  Fig.  A4. 

In  this  case  the  higher  order  coefficients  A^  and  A,-  are  set 
equal  to  zero.  When  all  the  coefficients  in  Eq.  A(13)  are 
known,  the  slope  factor  m,  for  any  given  (T  and  can  be 
calculated  from* 

m = ( a.  q + a^tf  + a2^  + a^<f3  + a^tf  + a^q”3)  A (14) 

+ (a0X  + axV  + a2\r2  + a^1^3  + a^1**4  + , , ,+ 

:;+(aQ5  + a13<y  + ag-^2  + 3^<p  + a43<j4  + a^5^5)^5 

This  procedure  resulted  in  a single  equation  (Eq.  A(14))for 
m in  which  the  coefficients  were  analytic  functions  of 
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The  coefficients  are  given  in  Tables  A1  and  A2.  Due  to  the 
well-behaved  nature  of  the  curves,  the  slope  factor  m could 
be  calculated  to  well  within  the  accuracy  of  the  original 
data  contained  in  Ref,  1, 

All  data  necessary  to  predict  performance  for  circular- 
arc  blading  was  reduced  to  polynomials  using  a similar 
procedure.  The  polynomials  were  included  in  a computer 
program  to  calculate  blade  element  performance  for  specified 
input  geometry  and  flow  conditions, 

A4,  COMPUTER  PROGRAM 

A*».l  Description.  Program  "NASA36"  was  written  for  the 
Hewlett-Packard  9830A  calculator  to  calculate  the  incidence 
and  deviation  angles  for  a rotor  with  meridional  inlet  ve- 
locity, The  program  follows  the  method  of  Ref,  1 through 
the  equations  given  in  Section  A2,  Where  a value  from  a 
Figure  in  Ref,  1 is  required,  it  is  calculated  using  a poly- 
nomial approximation  to  the  data  ss  described  in  Section  A3. 

The  program  is  self-explanatory  through  the  listing 
given  in  Table  A3.  The  symbols  are  defined  in  Table  A^, 

The  program  performs  the  blade  element  calculations  for  a 
specified  compressor  operating  condition.  The  inputs  are 
referred  flow  rate  and  speed,  and  temperatures  upstream  and 
downstream  of  the  rotor.  Through  the  method  given  in  Appen- 
dix C,  these  inputs  define  the  inlet  velocity  as  a function 
of  the  radial  position.  It  should  be  noted  therefore  that 
the  present  version  of  the  program  includes  a description 


of  the  transonic  compressor  inlet  flow-field,  while  it 
allows  changes  in  blading  geometry  to  be  input, 

A4.2  Operation.  Detailed  information  on  the  operating 
procedures  for  the  HP-9830  calculator,  mass  memory  and  as- 
sociated equipment  are  found  in  Ref.  9 and  Ref.  10.  The 
following  is  a step-by-step  operating  procedure  for  the 
program  called  "NASA36". 

1.  After  loading  "NASA36"  into  the  Hewlett-Packard 
9830A  calculator  press  "Run"  and  "Execute". 

2.  The  calculator  will  then  display  a series  of  inputs 
required  to  determine  the  geometry  of  the  blading.  These 
are  shown  in  Fig.  3 and  arei 

a)  Stagger  angle 

b)  Chord  with  sharp  leading  and  trailing  edges 

c)  Tip  radius 

d)  Hub  radius 

e)  Diameter  of  cone  through  stacking  radius 

f)  Maximum  Profile  Thickness 

g)  Radius  of  Profile  Contour 

h)  Blade  camber  angle 

i)  Profile  chord 

j)  Angle  of  cone  thru  stacking  radius 

3.  The  operating  parameters  of  the  compressor  rotor 
are  now  required.  These  arei 

a)  Revolutions  per  minute 

b)  Total  temperature  upstream  of  rotor 

c)  Temperature  downstream  of  rotor 


Coefficients  for  Eq.  A(ll) 


Table  A2 
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FIG.  A2  TOTAL-PRESSURE -LOSS  COEFFICIENT  VS.  INCIDENCE  ANGLE 


(IBIJBMI  tfftf-iTVfl 


SOLIDITY 

FIG.  A4  VARIATION  OF  COEFFICIENTS  OF  POLYNOMIAL  CURVES  IN  FIG.  A3 


APPENDIX  Bi  CALIBRATION  AND  CORRECTION  OF  PNEUMATIC  PROBES 


Bl.  FREE-JET  APPARATUS 

The  free-jet  apparatus  used  to  calibrate  probes  is 
shown  in  Fig.  Bl.  The  device  consisted  of  a 8 inch  diameter 
pipe  with  the  air  supplied  by  the  laboratory's  central 
Allis-Chalmers  Compressor.  A 4,2  inch  converging  nozzle 
was  used  to  increase  velocity  and  attain  Mach  numbers  up 
to  about  .9.  To  determine  if  the  exiting  velocity  was  uni- 
form and  undisturbed  by  the  probe  mounting  arrangement, 

surveys  were  conducted  with  a Prandtl  probe  over  a range 

X 

of  5 Mach  numbers  at. two  axial  locations.  Surveys  were 
also  conducted  with  the  same  probe  mounted  and  traversed 
from  opposite  sides  of  the  jet.  The  results  shown  in 

V 

Fig.  B2  indicated  that  the  velocity  profiles  were  uniform 
over  the  core  of  the  jet  extending  to  a diameter  of  3.8 
inches.  As  suggested  in  Ref.  11,  the  non-dimensional  velo- 
city (X)  appearing  in  Fig.  32  is  defined  by 

x = vAt  B{  1 ) 

where  v^5/^t  = . 62^T^  when  T^  is  in  °R.  Following 

Ref.  11  further,  the  relationships  for  stagnation  pressure, 
total  temperature,  and  Mach  number  are,  in  terms  of  X given 

by 

-J-  = 1 - X2  B(2) 

At 

-f-  = (1  - X2)  V(*-l)  B(3) 

rt 
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1 


X 

'1  - x< 


B(4) 


The  values  of  X in  Fig.  B2  were  calculated  using  the  meas- 
urement of  stagnation  pressure  and  the  assumption  of  atmos- 
pheric static  pressure  in  Eq.  B(3). 


B2.  CALIBRATION  OF  THE  COMBINATION  PROBE 

Calibration  of  the  pneumatic  probe  reported  in  Ref,  2 
can  be  represented  by  defining  two  functions  of  the  non- 
dimensional  velocity,  X,  and  the  pitch  angle, 0.  If 
and  ?2/3  are  indicated  stagnation  and  indicated  static  pres- 
sure of  the  probe  respectively,  then  define 


F©  = Fe  (X,0)  = 


P1  " P2/ 


B(5) 


P,  - P. 


Fpj  = Fj*  (X,0)  = — : 1-p- 


B(6) 


to  be  calibration  functions  when  the  probe  is  rotated  to 
balance  P2  = P^.  P^  is  the  fourth,  off  axis  pressure  meas- 
ured by  the  probe. 

The  complete  calibration  procedure  is  given  in  detail 
in  Ref.  3.  The  probe  was  recalibrated  to  cohfirm  the  re- 
sults given  in  Ref.  2 before  investigating  corrections  for 
boundary  effects.  A range  of  Mach  numbers  (X  = .13  to  .4) 
and  pitch  angles  (-20°  to  20°)  were  used.  With  true  static 
pressure  taken  to  be  atmospheric  and  true  impact  pressure 
that  measured  by  a reference  Prandtl  probe,  the  value  of 
X was  calculated  from  Eq.  B(3).  As  in  Ref.  2,  polynomial 


8: 


representations  were  found  for  the  calibration  equations 
written  in  the  form 

X = X(F0,0)  B(7) 

and 

F*  = Pf  (X,0)  B(8) 

where  F0  , and  are  the  ratios  of  probe  pressure  measure- 
ments defined  in  Eq.  B(5)  and  Eq.  B(6).  It  is  noted  that 
Eq.  B(6)  departs  from  the  definition  used  by  Dodge  in  Ref. 
2.  The  present  method  was  found  to  give  more  accurate  re- 
sults. The  coefficients  in  Eq.  B(7)  and  Eq.  B(8)  derived 
from  the  calibration  test  results  are  given  in  Table  Bl. 

As  a test  of  the  analytical  representation  of  the 
calibration,  the  calibration  data  were  used  to  calculate 
values  of  X and  0 to  compare  with  the  actual  values.  Table 
B2  shows  the  results  of  this  comparison  for  the  velocity. 
The  error  in  X is  generally  less  than  1 fo  and  only  reaches 
higher  values  at  high  positive  or  negative  pitch  angles 
( 16°  or  -8°)  or  at  higher  Mach  numbers  (.88).  At  this 
point  the  error  increases  to  about  2 %. 

Fig.  B3  shows  the  error  in  the  pitch  angle  in  the  • 
range  of  pitch  angle  from  -4°  to  +12°.  The  error  is  seen 
to  be  tl£°  and  is  not  a simple  function  of  either  pitch 
angle  or  Mach  number  (or  X).  There  appears  to  be  a limit 
to  the  capability  of  resolving  0.  The  error  was  larger  at 
high  negative  and  positive  values  of  pitch  angles  but  these 
were  outside  the  intended  range  of  use. 


B3.  BOUNDARY  EFFECT  ON  COMBINATION  PROBE  MEASUREMENTS 


In  application,  the  probe  is  not  in  plane  flow  away 
from  all  boundaries.  The  indicated  velocity  and  pitch  angle 
will  depend  upon  several  factors  including  boundary  layer 
and  stem  immersion  effects.  The  method  of  approach  assumed 
here  was  that  the  probe  always  obeys  the  calibrations  es- 
tablished and  that  errors  are  induced  due  to  the  perturbation 
of  the  probe  itself  on  the  flow  near  the  boundary.  Thus, 
in  any  flow  the  probe  measures  a velocity,  X,  and  a pitch 
angle,  (p , which  are  the  sum  of  the  actual  and  an  "induced" 
component. 

Three  surveys  were  conducted  in  a 10  inch  pipe  up  to 
the  maximum  attainable  Mach  number  of  .55*  The  results  of 
calculating  velocity  applying  the  free  jet  calibration  given 
above  are  shown  in  Fig.  B4. 

Introducing  subscripts  a(actual),  i( induced),  and 
m<measured),  the  velocity,  X,  may  be  represented  by 

= = B<9) 

where  the  fractional  'error*  in  X,  £x,  which  is  induced  by 
the  presence  of  the  wall,  is  given  by 

X - Xa 

= b(  10) 

X is  determined  from  the  probe  calibration  in  the  usual 
m 

way.  X is  the  required,  correct  velocity,  and  £ is  the 

CL  X 

error  to  be  established  by  calibration  tests. 


MB! 


r 


The  measured  velocity  is  corrected  using 

x --3s_ 

a ' 1 +£x 


B ( 1 1 ) 


when  the  error  £ is  known. 

In  the  same  manner  actual  pitch  angle  may  be  repre- 
sented by 

$a  = 0 - B(12) 

In  the  10  inch  pipe  tests,  the  actual  velocity  was 
calculated  from  the  probe  impact  pressure  and  pipe  wall 
static  pressure  using  Eq.  B(3)  and  the  actual  pitch  angle 
was  zero.  Thus  the  "induced"  error  in  X was  calculated 
using  Eq.  B(10)  and  the  "induced"  error  in  $ was  calculated 
using  Eq.  B(12).  The  two  errors,  and  , were  then  ex- 
amined as  a function  of  distance  from  the  pipe  walls  ex- 
pressed as  a number  of  probe  diameter?; (the  characteristic 
dimension  of  the  probe  was  taken  to  be  0.1  inches).  It 
was  found  that  the  errors  became  negligible  at  about  9.5 
probe  diameters  from  the  boundary.  In  order  to  derive  an 
expression  of  a suitable  form,  the  errors  were  plotted  as 
a function  of  (9.5-D),  where  D is  the  displacement  as  a 
number  of  probe  diameters.  The  velocity  error  was  then 
approximated  by  a polynomial  of  the  form 


€x  = A0  + Ax  (9.5  - D)  + A2  (9.5  - D)2  B(13) 

where  D is  measured  from  the  wall  through  which  the  probe 
was  inserted.  The  result  is  shown  in  Fig.  B5.  The  error 
was  less  than  2#  at  distances  greater  the  0.25  inches  from 
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the  wall.  In  order  to  extrapolate  the  correction  to  larger 
values  of  velocity,  the  error  was  represented  as  a function 
of  X,  The  coefficients  in  Eq.  B(13)  were  plotted  as  a func- 
tion of  X,  and  the  result  is  shown  in  Fig.  B6.  AQ,  and 
A2  were  seen  to  be  almost  linear  functions  of  X and  were 
therefore  approximated  by 

A0  * B0  + B1X,95  + B2X.952 

Al  = co  + C1X,95  + C2X.952 

A2  * D0  + D1X,95  + D2X,952 

where  X ^ was  the  value  of  X at  D = 9.5«  Values  of 
X(corrected)  minus  X( uncorrected)  verses  immersjon  depth  are 
shown  in  Fig.  B7.  Because  of  the  aerodynamic  design  of  the 
probe,  the  corrections  necessary  for  boundary  effects  are 
small. 

Results  for  the  induced  pitch  angle  error  are  shown  in 
Fig.  38.  Since  these  errors  were  within  the  resolution 
capability  of  the  probe  (l£°),  no  correction  was  made. 

B4.  COMBINATION  PROBE  MEASUREMENTS  IN  STEADY  SWIRLING  FLOW 
Additional  errors  may  be  present  when  a pneumatic  probe 
is  immersed  in  a swirling  flow.  Fig.  B9  shows  the  apparatus 
which  was  constructed  to  investigate  this  possibility.  The 
apparatus  consisted  of  an  annulus  with  a channel  height  of 
l£  inches  (10  inches  in  outer  diameter).  A total  of  40 
static  wall  readings  were  available  and  a provision  was 
made  for  probe  surveys  at  different  axial  stations.  The 


swirl  velocity  was  generated  by  a single  stage  of  stator 
blading  at  the  entrance  to  the  annulus.  Initial  measure- 
ments suggested  the  possibility  of  boundary  layer  separation 
and  reverse  flow  near  the  hub.  In  an  effort  to  correct  this 
problem,  a plate  was  constructed  to  decrease  the  channel 
height  at  the  exit,  isolating  the  pressure  at  the  hub  wall 
from  the  surrounding  atmosphere,  and  blocking  backflow. 

Results  are  given  in  Pig.  BIO  which  shows  the  static 
pressure  derived  from  the  probe  calibration.  Agreement 
is  very  close  near  the  hub  with  a moderate  departure  near 
the  tip  at  the  highest  Mach  number. 

B 5.  CALIBRATION  OF  THE  WEDGE  PROBE 

A type  of  probe  which  is  reported  to  give  a good  indi- 
cation of  static  pressure  over  a wide  range  of  Mach  number 
is  the  wedge  static  probe  illustrated  in  Fig.  5.  The  probe 
has  a hatchet-shaped  head  with  small  pressure  taps  in  each 
face  of  the  wedge.  In  principle,  if  the  leading  edge  is 
sharp  and  the  wedge  angle  is  small,  the  flow  past  the 
sensing  orifice  will  be  almost  unperturbed,  and  the  meas- 
ured pressure  will  be  close  to  the  static  pressure  in  the 
undisturbed  flow. 

The  static  pressure  may  thus  be  determined  by 

Ps  - < P2/3  - PA  > + PA  B<U> 

where  it  is  assumed  the  probe  is  aligned  in  the  flow  with 
and  equal,  and  PA  is  the  reference  atmospheric 
pressure. 
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The  wedge  probe  shown  in  Fig.  5 was  traversed  in  both 


the  free  jet  and  the  swirling  annular  flows  through  the 
same  range  of  Mach  numbers  reported  for  the  combination 
probe.  However,  the  characteristics  of  the  probe  in 
swirling  flow  were  not  fully  determined.  It  was  found  in 
analyzing  the  results  of  the  tests  in  the  swirl  annulus 
that  the  radial  gradient  in  the  flow  angle  gave  rise  to 
probe  readings  which  were  in  error  when  the  probe  was  ro- 
tated to  balance.  This  was  explained  by  the  difference  in 
the  radial  position  of  the  sensors  on  the  two  faces  of  the 
wedge.  It  was  therefore  necessary  to  use  the  combination 
probe  to  sense  flow  direction  in  the  compressor  measure- 
ments and  to  set  the  wedge  to  the  measured  angle. 
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Table  B2 

Error  Analysis  of  Limiting 
Velocity  From  Combination  Probe  Calibration 


4> 

foE 

X 

X1 

%E 

X 

X1 

-8 

1.240859 

.131396 

.133026 

1.294396 

.164529 

• 166655 

-4 

. .179254 

.131396 

.131631 

0 

.278352 

.131396 

.131761 

.500537 

.164718 

.165542 

4 

-1.363128 

.131396 

.129604 

1.408336 

.164718 

.01670; 

8 

-1.530939 

.131396 

.129384 

.203496 

.164718 

.16505; 

12 

-0.680858 

.131396 

.130501 

.651461 

.164718 

.165791 

-8 

1 . 000402 

.19477 

.196718 

.233878 

.220488 

.22100: 

-4 

-0.117261 

.194314 

.194086 

.152853 

.220744 

.221081 

0 

.169147 

.104618 

.194947 

-.114982 

.220616 

.220362 

4 

-0.214012 

.194161 

.193745 

-.092819 

.221253 

.22104? 

8 

-0.358790 

.194466 

.193768 

-1.029916 

.220871 

.218595 

12 

.376707 

.19477 

.195503 

-.646592 

.220616 

,21918s 

-8 

-1.239331 

.246923 

.243862 

.237152 

.287457 

. 288135 

-4 

-0.034597 

.247245 

.247159 

1.076403 

.287457 

.290551 

0 

-0.146786 

.247138 

.246775 

1.024665 

.28754 

.290486 

4 

-0.515128 

.24703 

.24 5757- 

.471649 

.287872 

.28922$ 

8 

-1.285034 

.24703 

.243 855 

-1.434447 

.287872 

.283742 

12 

-1.132503 

.246923 

.244126 

.699168 

.287872 

.2898814 

-8 

4.327326 

.330596 

.344901 

.274718 

.342723 

.3436614 

-4 

-3.051296 

.330723 

.340814 

-1.857185 

.342664 

.33630C 

0 

1.741821 

.329579 

.335319 

-1.009610 

.34284 

.339375 

4 

2,010646 

.33034 3 

.336985 

-.969406 

.342488 

.33916? 

8 

2.250762 

.330089 

.337518 

-1.002690 

.342488 

.33905; 

12 

-0.449206 

.330089 

.328606 

-1.115949 

.342547 

.3387214 

0 = Pitch  Angle 

X = Actual  Value 

X^  = Computed  Value 

%E  = Error 

SIDE  VIEW 


FIG.  B1  FREE-JET  APPARATUS  USED  TO  CALIBRATE  PROBES 
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APPENDIX  C » INLET  FLOW  FIELD  DETERMINATION 


Cl.  METHOD  OP  APPROACH 

An  important  consideration  in  a transonic  compr  ssor 
installation  is  the  velocity  profile  at  the  face  of  the 
rotor.  The  profile  depends  on  the  inlet  shape  and  the  mass 
flow  rate.  The  following  analysis  is  based  on  the  assump- 
tion that  a close  approximation  to  the  true  velocity  may 
be  calculated  if  the  referred  flow  rate  is  known. 

Seven  surveys  were  conducted  by  Anderson  (Ref.  3)  in 
which  accurate  calculation  of  the  inlet  velocity  profiles 
was  accomplished  using  a calibrated  United  Sensor  Probe 
over  a range  of  flow  rates  tabulated  in  Table  Cl. 

Reference  11  defines  a "total  flow  function,"  , given 
by 

X 3 -IT-  = FT"  = F(x>  = XU  " x2)1/fir"1)  C(l) 

The  flow  function  provides  the  means  of  replacing  density 
in  the  equation  with  the  non-dimensional  velocity  (X)  for 
internal  compressible  flows.  For  each  survey  conducted  in 
Ref.  3»1Ewas  calculated  at  each  point  and  plotted  versus 
the  normalized  displacement  (hub  to  tip).  The  results  are 
shown  in  Fig.  Cl.  Each  curve  was  then  fitted  with  a suit- 
able polynomial  of  the  form 

jj[  * aQ  ♦ atH  + a2H2  + a3H3  + ...  C(2) 

where  H * normalized  displacement  (hub  to  tip)  and  a1 , a2,.. 
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are  determined  by  a polynomial  regression  curve  fitting 
routine.  The  value  of  the  total  flow  function  at  the  cen- 
terline was  evaluated  for  each  case.  Values  of 

^ici,  for  each  survey  were  then  calculated  and  plotted 
against  H.  The  resulting  curves  shown  i'n  Pig.  C2  were  then 
approximated  by  a family  of  polynomials  having  the  form 

5 =5/£Cl  * Aq  + AtH  + A2H2....  C(3) 

\ / 
which  is  then  called  the  normalized  flow  function. 

The  referred  flow  rate  is  defined  as 


where 


* tbYe 

" S’ 


C(4) 

C(5) 


and 


and  O = mass  flow  rate, 

A plot  of  JCL  versus  sfc  wn  in  Fig.  C3  yielded  the 
following  linear  relationship! 


V- 


t0 


C ( 6 ) 


Thus, 


IcL»  bQ  + b^*  C(7) 

Eq.  C(3)  and  Eq.  C ( 7 ) provide  an  analytical  rep- 


resentation of  the  inlet  flow  profile  in  terms  of  the 
referred  flow  rate  only. 


C 2.  APPLICATION 

For  each  run  of  the  transonic  compressor,^  is  calcu- 
lated and’jJT  is  then  evaluated  from  Eq.  C(7).  can  then 
CL 

be  evaluated  from  Eq.  C(3). 

X is  then  evaluated  as  the  subsonic  solution  of  the 
flow  function  given  by  Eq.  C(l)  using  Newton's  method  as 
follows i 

1)  Guess  X = XQ 

2)  Calculate 

3>2  = X d - x2)1/(Y”1)  C(8) 

lx"  = ^2  X ‘ (y-i)  ( 1_XJC2)  c(9) 

2 

X = x + 0(10) 

dX 

3)  Now  test  for  convergence  of  3E . 

If  \12-1  ^|<. 00001  then  the  calculation  is  complete. 
If  this  requirement  is  not  satisfied  then  the  new  in 
Eq,  C(8)  is  calculated  using  the  new  X in  Eq.  C(10), 

4)  Mach  number  is  calculated  using  Eq.  B(4), 


Table  Cl 

Flow  Rates  Used  for  Inlet  Velocity  Determination 


Referred 
Flow  Rate 

Centerline  Value 
of  Flow  Function 

8.95 

.09589 

9.36 

.10160? 

10.89 

.11693 

11. lS 

.12039 

11.42 

.122025 

13.83 

.151239 

14.30 

.156789 

NDN-DJMEN51DNHL  DISTANCE  (HUB  TO  TIP) 

PIG.  Cl  TOTAL  FLOW  FUNCTION  VS.  NORMALIZED  DISPLACEMENT 
UPSTREAM  OF  THE  ROTOR 


NDRHflUZED  DI5PUKENENT  (HUB  TO  TIP) 

FIG.  C2  NORMALIZED  FLOW  FUNCTION  VS.  NORMALIZED  DISPLACEMENT 
UPSTREAM  OF  THE  ROTOR 


APPENDIX  Di  ROTOR  LOSSES 


Dl.  NASA  LOSS  COEFFICIENT 

Losses  incurred  in  flows  in  turbomachinery  may  be  de- 
scribed with  a varity  of  methods  depending  upon  the  measure- 
ments available  and  the  specific  application.  Ref.  1 
describes  a loss  coefficient  defined  by 


o'-  ■■■■  »<*> 

*1  " n 

where 

^ = mass-averaged  defect  in 

relative  total  pressure 
divided  by  the  pressure 
equivalent  of  the  inlet 
velocity  head 


P.  = relative  total  pressure 
upstream 

P ' = relative  total  pressure 
^ downstream 

P^.  = static  pressure  upstream 

9 

(P-).,  = ideal  relative  total 
d 1Q  pressure  downstream 

From  the  measurements  obtained  and  a knowledge  of  the 

inlet  velocity  profile  the  parameters  P^,  P2,  P.^  ^t2*  *1* 

Y m m —I 

2*  ti*  t2  are  known.  A representation  of  W in  terms 
of  measured  stagnation  conditions  and  velocities  (X  or  M) 
is  desired. 
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Dl.l  Upstream  of  the  Rotor.  The  following  relations 
may  be  writteni 


<Vid  * p 


tl 


*/*-!) 


y2  1 */*-!) 

XU2J 


D(2) 


D(3) 


Tti  = 1 " X1  + + XU2  “ XU1 

where  X is  a velocity  non-dimens ionalized  as  in  Eq.  B ( 1 ) , 
subscript  denotes  the  relative  velocity,  and  subscript  W 
denotes  the  wheel  speed.  The  other  notation  is  given  in 
the  T-S  diagram  for  the  flow  through  the  rotor  shown  in  the 
above  sketch. 

Using  the  above  equations,  for  an  "impulse"  machine 
for  which  there  is  a meridional  entrance  flow,  the  following 
relations  result i 


■tl 


= 1 + X 


U2 


2 n */(*-!) 


and 


fi  - ?ti  < 1 + x5i  > 


P1  ’ Pt!  ( 1 - XP  ) 


D(5) 

D(6) 

D(7) 


D1.2  Downstream  of  the  Rotor.  Sufficient  measurements 

t 

exist  at  the  machine  station  2 to  derive  ^2  using  Tt2#  T^2, 
X2,<*2,  and  U2  where  X2  is  the  measured  downstream  velocity, 
non-dimensionalized  using  the  downstream  limiting  velocity 
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(V^g)  in  Eq.  B(l).  Alternatively,  since  Tg=  constant  for 
the  flow  through  the  rotor  (see  diagram),  then  one  of  the 
separate  measurements  is  not  needed.  Two  methods  were 
therefore  used.  Method  1 used  all  measurements  at  station 
2.  Method  2 did  not  use  the  temperature  rise  at  station  2, 
and  thus  served  as  a check  on  the  accuracy  of  the  temper- 
ature measurement. 

(i)  Method  I (all  measurements  at  station  2) 


p.  - P fIlo*/(*- 1) 

- rt2  ) 


t2 


D(8) 


where  * — is  given  by 
1 12 


r2  ,Tti, 

” * ' AJJZ  ~ *"(ii) * * v2  ~UaV  Tt2  OJ-“^2 

Using  Eq.  D(6),  D(7)»  D(8),  and  D(9)  in  Eq.  D(l),  there 
results 


= 1 - (nr*)  - 2X,  xttV!L  Sin<x, 


•tl 

v2  ~UaV  Tt2 


D(9) 


CO  = 


fl  + X2  -ll  - x21 

l uij  L AiJ 


*U2  2X2Xu2sir^>  yiY  ^tj. , 


[i+Xy y/(V_1)  - 


D(  10 ) 


(ii)  Method  II  (Temperature  measurement  at  station  2 

~ not  needed) 

Temperature  measurement  downstream  may  be  found  using 

the  following  relation. 


XE  _ 
*t2 

Define 


- (-rfj-J/dtaAti)  - d+x|2)/(Tt2/rtl) 


■n  r Tt2 

Y = 


X2Xu2  sinOCg  +y!72X U2  sin  z)c  + .1 


D(  11 ) 


D(12) 
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Using  Eq.  D(6),  D(7),  D(ll),  and  D(12)  in  Eq.  D(l),  the 
loss  coefficient  is  given  "by 


,n.  KJ  V(v-l)-  <£f) 

" [l  ♦ x§JV/'K-1)  -[l  - xf] 

D2. . DIFFUSION  FACTOR 


D(13) 


Ref.  1 reports  results  of  correlating  compressor  losses 
in  terms  of  the  "diffusion  factor".  This  approach  results 
from  the  argument  that  the  loss  in  stagnation  pressure  de- 
pends largely  on  the  suction-surface  adverse  pressure  gra- 
dient, of  which  the  diffusion  factor  can  be  shown  to  be  an 
approximate  measure.  The  diffusion  factor  of  Ref.  1 is 
defined  by 


D = 1 


W2  RiwUi  - R2WU2 

W1  & (Rx  + R2)W1 


D(14) 


In  this  form,  changes  in  axial  velocity  and  radius  through 
the  blade  row  are  accounted  for. 

It  is  concluded  in  Ref.  1 that  the  total-pressure-loss 
parameter  (fl'c2^  -2)  can  be  correlated  in  terms  of  the 
above  diffusion  factor.  Each  family  of  blade  shapes  has  a 
somewhat  different  correlation  curve,  however,  this  pre- 
sents no  real  problem  as  the  total  number  of  shapes  is 
limited.  The  correlation  for  several  types  of  blading  is 
shown  in  Fig.  203  of  Ref,  1. 
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